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Abstract 



It is pointed out that "the colorless objects(s)" in diffractive lepton-nucleon scat- 
tering in the small-xs region can be probed by examing the final state particles 
in such processes: (A) Perform multiplicity and/or energy factorial moment analy- 
ses, and check especially the xp- and Q 2 -dependence of the intermittency- indices. 
(B) Examine the diffractive structure function F^ 3 \(3 , Q 2 ; x p) as function of xp 
for different (3 and Q 2 . It is shown that such analyses and measurements can yield 
useful information on the hadronization process and on the intrinsic properties of 
the colorless objects. In particular, the observation of power-law distributions of 
kP^(/3, Q 2 ;xp) with respect to ip and with respect to xp — independent of (3 and 
Q 2 -values may be considered as strong indication for the existence of self-organized 
criticality in gluon-cluster formation processes. 







The observation of large rapidity gap events in deep-inelastic lepton-nucleon scattering 
experiments performed in the small xp region§ shows that the virtual photons in such pro- 
cesses may encounter colorless objects originating from the nucleon. What are such objects? 
While the existing datJ can be reasonably well reproduced by Regge-pole mo deli, it is not 
known whether/how the properties and the formation process(es) of such objects can be 
directly probed — - in a model-independent manner. Furthermore, deep-inelastic lepton- 
nucleon scattering experiments and empirical analyses in this kinematical regionB! show 
that the gluon-density in the nucleon is much higher than those for quarks/antiquarks, and 
it is increasing with decreasing x^i. Are the striking phenomena observed in the small xp 
region related to one another? In particular, can the above-mentioned colorless objects be 
virtual states formed by interacting " soft" gluons which exist in abundance in this kinemat- 
ical region?! Can the observed final state particles tell us something about the space-time 
properties of the colorless object? Is it possible to extract such information directly from 
experiments? 

In the present note we show that these questions can be answered in the affirmative. 
Because of the complexity of such process, we propose to use concepts and methods of 
statistical physics and study distributions, average values, and fluctuations of the relevant 
observables which describe the final-state particles in the appropriate kinematical regions. 
To be more precise, we suggest the following: (A) Perform multiplicity (and/or energy) 
factorial moment analyses to obtain the indices of intermittency for different Q 2 - and re- 
values. (B) Examine the diffractive structure function F% (ft, Q 2 ; xp) as function of xp for 
different values of (3 and Q 2 . 

In the analyses mentioned in (A), we do the following: (i) Consider event-by-event, at 
given values of Q 2 and xp, the data for diffractive lepton-nucleon scattering in which the 
momentum-fraction (often denoted by xp) of the exchanged colorless object is known. The 
latter quantity can be measured by detecting the final state of the outgoing proton (or 
nucleon-resonances which have similar quantum numbers as the proton) or by determining 
it approximately by measuring Q 2 and M 2 in such processes, (ii) Focus attention on the 
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multiplicities of the charged hadrons produced in the collision between the virtual photon 
and the colorless object. That is, look at the multiplicities of the hadronic system which 
yields the " M x -spectrum" , and use them to calculate the normalized factorial moments 
(FM) F q as functions of the resolution of the phase space. To be more precise, consider the 
multiplicities of the hadrons in a given phase space interval characterized by A77, Ap±, A<f>, 
(pseudorapidity, transverse momentum and azimuthal angle of the observed particle). Divide 
the above-mentioned phase space into M subintervals and calculate! 

p 1 ^ {n m {n m -!)■■■ (n m - q + 1)) m 

where n m is the multiplicity in the mth subinterval, and (• ■ ■) means taking the average 
over events. Use the obtained results to calculate the intermittency index ip q , from the F q 
vs M plots, (iii) Repeat the procedure mentioned in i and ii for a fixed Xp, a fixed Q 2 , 
but for different xp values, (iv) Repeat the procedure mentioned in i and ii for a fixed (or 
integrated) xp, a fixed Xp, but for different Q 2 values. 

The proposal mentioned in (A) is based on the following reasonings: 
Factorial moment analyses can be used to extract information on dynamical processes in 
a model-independent way. In particular, if dynamical fluctuations have a typical size (e.g. 
5t]q in one-dimensional pseudorapidity phase space due to resonance decays which have a 
width Srjo = 1) the factorial moments rise with decreasing bin-size in pseudorapidity-interval 
Ar]/M as long as Arj/M > Srjo and saturate to some constant value for Arj/M < 5r]o. But, 
if self-similar fluctuations exist at all scales of Aij/M, the factorial moments are expected to 
follow a power-law F q ~ M^ q in the limit M — > 00 which manifests itself in the corresponding 
In F q vs In M plot as a linear rise of the slope. Such a behaviour is known!'! as intermittency, 
and its strength is characterized by the slope parameter <p q which is also called the anomalous 
scaling index. 

The xs-dependence of such F q vs M plots, and in particular that of (p q (Q 2 ,xp]Xp) for 
fixed Q 2 (and xp), plays a distinguished role in such studies. This is because, viewed from a 
fast moving reference frame, for example the lepton-nucleon center-of-mass frame, where the 



nucleoli's momentum P is large in high-energy collisions, we see in a diffractive scattering 
eventB the following: The virtual photon 7* originating from the incident lepton is absorbed 
by the nucleon. Precisely speaking, its energy- momentum q = (q°, q) is absorbed by a virtual 
colorless object of the nucleon. The time interval T- in t in which the absorption process takes 
place (it is known as the lepton- nucleon interaction/collision time) can be estimated by 
making use of the uncertainty principle. In fact, by calculating l/q° in this reference frame 
we obtain: 

4IPI x B 
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Q 2 \-x B 

This means, for given P and Q 2 = — q 2 , r^t is directly proportional to xb for xb « 1- 
In other words, Xb is a measure of the time-interval in which the absorption 0/7* by the 
space-like virtual colorless object takes place. Hence, by studying the x^-dependence of the 
intermittency index ip q (xB, Q 2 ', xp), we are not merely probing the statistical and dynamical 
fluctuations of the collision process between 7* and the colorless object which we hereafter 
call Cq. Since this process of hadronization of the virtual colorless object Cq is initiated by 
the interaction with 7*, we are also examing whether (if yes, how) the hadronization process 
changes with the interaction time r int . This question is of considerable interest, because a 
virtual photon 7* can (logically) only be absorbed by virtual systems (cjj's) whose lifetimes 
(t*'s) are longer than r int (i.e. r int < r c ). That is, the average lifetime (r c ) of the cjj's, which 
can absorb a 7* associated with interaction-time r int , is a function of the number of Cq's which 
satisfy the condition r int < r c . Hence, from the ^-dependence of the scaling behaviour of 
FM's, in particular from that of the corresponding <y9 g 's, we can find out whether/how the 
process of hadronization of a Cq depends on its average lifetime (r c ) of the Cq's. This means, 
by measuring the above-mentioned ^-dependence we can in principle obtain information 
on the number-distribution of such colorless objects as a function of their lifetimes. It should 
be mentioned at this place however that, in many cases (an example will be given below) 
multiplicities depend first of all on the invariant mass of the produced particles. Hence, 
special care has to be taken to conjecture such associations. See in this connection, also the 



method proposed in (B) below. 

The Q 2 -dependence of the scaling behaviour of the FM's, in particular that of 
ip q (Q 2 ,XB',xp), is also of considerable interest. This is because, in photon-proton scat- 
tering experiments, not only those with real (Q 2 = 0) photons but also those with space-like 
(Q 2 > 0) photons where Q 2 is not too large (< lGeV 2 /c 2 , say) have very much in common 
with hadron-hadron collisions. Having in mind that the index of intermittency for hadron- 
hadron scattering is smaller than that for electron-positron annihilation processes!, we are 
led to the following questions: Do we expect to see a stronger Q 2 -dependence when we in- 
crease Q 2 from zero to 10 or 100 GeV 2 /c 2 , say? Is this also a way to see whether space-like 
photons at large Q 2 "behave like hadrons" in such interactions? 

While waiting for data to perform the above-mentioned analyes, one may want to use 
some existing phenomenological models, for example JETSET0 to generate such " data" , and 
use them to obtain F q vs M plots, and to obtain the xb- as well as the Q 2 -dependence of 
the intermittency indices for fixed Q 2 and xb respectively. The results of such a calculation 
are predictions based on the assumption that JETSEli is applicable to diffractive lepton- 
nucleon scattering. These prediction can be tested when the data are available. By carrying 
out such model-calculations, we convinced ourselves that the proposed method works well. 
The details on this and on other model-calculations will be discussed elsewherei. 

Furthermore, having in mind that jets have been observed (see e.g. Ref.0 and the pa- 
pers cited there.) in diffractive electron-proton scattering processes, and experimentally 
calorimeters have been used to measure the collision events in general, and to study jet- 
events in particular, we are naturally led to ask: Can we, instead of using the distribution 
of multiplicities in phase space, also use the distributions of transverse-energies to probe the 
existence of dynamical fluctuations in such intermittency-analyses? Here, the transverse- 
energy E± is measured on an event-by-event bases with respect to the axis of the virtual 
photon. According to the most recent experimental knowledge (See e.g. the review paper 
given in Ref.l), we expect to see that the distributions of E± in phase space in such collision 
events are symmetric with respect to this axis, and symmetric with respect to the origion 



of the c.m.s. frame of the colliding objects — -namely 7* (the virtual photon) and Cq (the 
virtual colorless object). As a first step, we generalized in a straightforward manner the 
usual procedure^ by introducing an energy-unit e and write the "energy factorial moment" 
as (E±(E± -e)---[E ± -(q- l)e})/ (E ± ) q . It is clear that Eje can be considered as 
integers, provided that e is sufficiently small. Under this condition, it is clear that statistical 
fluctuations in F^ can be canceled out in the same way as that in F q defined in Eq.( |l]). 
But, this means, there is a dependence on an arbitrary parameter e, when we use F^l Is 
there a way to get rid of this kind of arbitrariness in the practice? 

In order to answer this question, let us consider a collision event, in which the total 
transverse energy is Ex (total), and the ratio between the arbitrarily chosen e and E± (total) 
is A i.e. A = s/E± (total). We generate events under the assumption that E±/e obeys the 
Bernoulli distribution for different sizes of subintervals, that is, when the above-mentioned 
number M changes. It is clear that the slope in the double logarithmic F^ vs M plot has 
to be flat. Next, we introduce the moment R q E > = (E\_) / (E±) Q and calculate this quantity 
for different choices of A. The corresponding R% vs M plots is shown in Fig.l. Here we see 
that i?2 can be considered a good approximation for F 2 , when e is of the order of 10~ 3 
of the total E± in the event under consideration. 

Let us recall that, by studying F q we are examing the anomalous scaling behaviour! of 
the probability-moments for multiplicies C q = 1 / M J2m=i(Pm) / (Pm) q ] by measuring R q E \ we 
are looking at that of the corresponding probability-moments for transverse energies. Since 
the R q E > vs M plots for transverse energies distributions are in general different from the 
F q vs M plots for multiplicities, comparisons between the two kinds of plots, together with 
the xb- and Q 2 -dependence of the corresponding intermittency indices, will raise a number 
of questions, in particular the following: Is the multiplicity distribution or the transverse- 
energy distribution more relevant for studying dynamical fluctuations? It seems that much 
work is still needed to answer such questions. 

As we have mentioned in the introduction, and have shown in the part (A) of this paper, 
our goal is to find out whether/how model-independent information on the colorless objects 
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can be obtained by examing the final state particles in diffractive lepton-nucleon scattering 
processes. For this purpose, we now consider the measurements mentioned in (B), and discuss 
the notion as well as the consequencies of "self-organized criticality" . We recall that, in a 
series of papers, P. Bak, C. Tang and K. Wiesenfeld@(hereafter refered as BTW) have shown 
that certain dissipative dynamical systems with extended degrees of freedoms can evolve 
towards a self-organized critical state which give rise to spatial and temporal power-law 
scaling behaviour. The spatial scaling leads to self-similar fractal structure. The temporal 
scaling manifests itself as flicker noise. Having in mind that the density of soft-gluons is high 
in the small xb region!, and that in an extended system (in space-time) such gluons may 
interact with one another in accordance with QCD, the facts found by BTW leads us to 
the question: Can systems of gluons of the target proton evolve into self-organized critical 
states? Knowing (at the present time) very little about the dynamics of soft-gluons on the 
one hand, and even less about the general dynamical origin of self-organized criticality on 
the other, we are not in a position to answer this question theoretically. But, with the 
help of our colleagues at HERA, there seems to be a chance to answer it experimentally. 
To be more specific, let us rephrase the question as follows: If the interacting soft-gluons 
can indeed form such critical states, where the interaction with an additional gluon can be 
considered as a local perturbation, what kind of signals do we expect to see experimentally? 

According to the known characteristics of self-organized states of BTW©, the structures 
of such states are barely stable; and a local perturbation of a critical state can grow over all 
length scales, leading to anything from a "shift" of a single unit to an avalanche, the size 
of which can be as large as the entire system. The distribution of the transported physical 
quantities, "the dissipated energy" , is a measurable quantity; and it has been shown by BTW 
(in two and three dimensional cases) that the distribution of this energy (also called "the 
size of cluster" by BTW), and the distribution of the time- interval in which the formation 
process of such a cluster takes place (called by BTW "the lifetime of the cluster") always 
obey power-law. 

In connection with the expected experimental signals, let us recall the following: Diffrac- 
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tive lepton-proton scattering can be envisaged as the collision between the virtual photon 
7* and the virtual colorless object Cq which carries a fraction xp of the energy of the in- 
cident proton (in, e.g. the above-mentioned electron-proton cm. frame). The "diffractive 
structure function" F® ( ' 3 \i3,Q 2 ]Xp) at given (3- and Q 2 - values can be considered as the 
probability for 7* (characterized by Q 2 and (3) to encounter the Cq. Hence, F®^ is simply 
the probability for the above-mentioned 7* to " which carries the energy fraction 

xp (in other words, "the xp-distribution of the Cq in the proton" if it is independent of /3 
and Q 2 ). This means, for given constant values of (3 and Q 2 , we may check whether the 
Cg's can be identified with BTW-clusters in the following way. Since in that case: (a) the 
interacting soft-gluons can evolve into a state with no characteristic time or length scale, 
(b) the gluon-clusters, Cq s, circ the results of perturbation of a self-organized critical state 
caused by soft-gluon interaction, and (c) "the distributions of the dissipative energy" (of 
BTW) are nothing else but the distributions of energies of the gluon-clusters, we should 
see that F^ 3 \f3, Q 2 ;xp) obeys a power-law Xp a with a ~ 1, for given constant values of (3 
and Q 2 . Furthermore, keeping in mind, (as we have already discussed in connection with 
proposal A) that 15 is a measure of the average lifetime of gluon-clusters; and that the 
variables Xb, (3 and xp are related to one another by xp = (3xp, we also should see that 
P^F^ 3 ' ((3,Q 2 ;xp) show a power-law behaviour with respect to xp (Note that (3 can be 
considered as the fraction of the 4-momenta of Cq carried by the struck charged constituent 
interacting with 7*). Let us now look at the experimental datJ and examine first the xp- 
dependence of F^ (/3, Q 2 ; xp) at given values of (3 and Q 2 . As we can see in Fig.2 ( Note 
that the corresponding plots for f3~ 1 F^ 3 ^ vs xp = f3xp have exactly the same form for 
given f3\), the existing datJ indeed show the expected characteristic spatial and temporal 
scaling behaviours. These results may be considered as strong indication that self-organized 
criticality plays a dominating role in the formation of gluon-clusters, the color-singlet ones of 
which are the objects we are dealing with in diffractive lepton-uncleon scattering processes. 
Further implications of this observation will be discussed in more detail elsewhere!. 
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Figure captions 

Fig.l, lni?f = \n({E'j_)/{E_ L } 2 ) as functions of InM = ln(A/<5), when the transverse energy 
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E± in units of e in subinterval 5 is stochastically produced according to Bernoulli distribu- 
tion. Here, A = e/Ej_ (total). 

Fig. 2, The dependence of F2 ((3,Q 2 \ xp) upon xp at given values of (3 and Q 2 . The data 
are taken from Refill (The data at lower Q 2 are omitted because of lack of space. They also 
fall on one straight line in such log-log plots consistent with the data at large Q 2 values). 
See text for further details. 
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